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>«.i^ ■ We consider a wireless node that randomly receives data from different sensor units. The arriving 

Cn ' data must be compressed, stored, and transmitted over a vifireless link, where both the compression 

and transmission operations consume power. Specifically, the controller must choose from one of 
multiple compression options every timeslot. Each option requires a different amount of power 
and has different compression ratio properties. Further, the wireless link has potentially time- 
^+ ' varying channels, and transmission rates depend on current channel states and transmission power 

fNl ' allocations. We design a dynamic algorithm for joint compression and transmission, and prove that 

it comes arbitrarily close to minimizing average power expenditure, with an explicit tradeoff in 
average delay. Our approach uses stochastic network optimization together with a concept of place 
C J ' holder bits to provide efficient energy-delay performance. The algorithm is simple to implement 

/^*s ' and does not require knowledge of probability distributions for packet arrivals or channel states. 

^^ Extensions that treat distortion constraints are also considered. 

' ■^ ' Categories and Subject Descriptors: ... [...]: ... 

C^ , General Terms: queueing analysis, stochastic network optimization, sensor networks, data fusion, 

distortion, fading channel 



00 ; 1. INTRODUCTION 

r~^, , We consider the problem of energy-aware data compression and transmission for 

ff^ ■ a wireless link that receives data from N different sensor units (Fig. 1). Time is 

slotted with normalized slot durations t € {0, 1,2,.. .}, and every timeslot the link 
(^ ' receives a packet from a random number of the sensors. We assume that packets 

00 . arriving on the same timeslot contain correlated data, and that this data can be 

^^ ' compressed using one of multiple compression options. However, the signal process- 

ing required for compression consumes a significant amount of energy, and more 
. ; ■ sophisticated compression algorithms are also more energy expensive. Further, the 

r> I data must be transmitted over a wireless channel with potentially varying channel 

j^ ' conditions, where the transmission rates available on the current timeslot depend 

on the current channel condition and the current transmission power allocation. 
The goal is to design a joint compression and transmission scheduling policy that 
minimizes time average power expenditure. 
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Fig. 1. Multiple sensors sending data to a single wireless link. 



This problem is important for modern sensor networks where correlated (and 
compressible) data flows over power limited nodes. Compressing the data can save 
power by reducing the amount of bits that need to be transmitted, provided that 
the transmission power saved is more than the power expended in the compression 
operation. It is important to understand the optimal balance between compres- 
sion power and transmission power. Work in [Barr and Asanovic 2003] considers 
this question for a wireless link with fixed transmission costs, and describes prac- 
tical compression issues and reports communication-to-computation energy ratios 
for popular algorithms. Work in [Sadler and Martonosi 2006] considers a simi- 
lar static situation where the wireless channel condition is the same for all time. 
There, it is shown experimentally that compression can lead to a significant power 
savings when data is transmitted over multiple hops. The proposed algorithm of 
[Sadler and Martonosi 2006] uses a fixed data compression scheme, an adaptation 
of the Lempel-Ziv- Welch (LZW) compression algorithm for sensor networks. Tech- 
niques for distributed compression using Slepian-Wolf coding theory are considered 
in [Pradhan et al. 2002] [Cristescu et al. 2005]. Models of spatial correlation be- 
tween data of different sensors are proposed in [Pattem et al. 2004] and used to 
construct and evaluate energy-efficient routing algorithms that compress data at 
each stage. 

The above prior work has concentrated on static environments where transmission 
power is directly proportional to the number of bits transmitted and/or traffic 
rates are fixed and known, so that compression and transmission strategies can 
be designed in advance. Here, we focus attention on a single link, but consider a 
stochastic environment where the amount of data received every slot is random, as 
is the current channel condition for wireless transmission. Further, the transmission 
rate is an arbitrary (possibly non-linear) function of transmission power. Optimal 
policies in this stochastic context are more complex, and more care is required to 
ensure transmissions are energy-efficient. 

In this paper, we design a dynamic compression and transmission scheduling 
algorithm and prove that the algorithm pushes total time average power arbitrarily 
close to optimal, with a corresponding tradeoff in average delay. We assume the 



algorithm has a table of expected compression ratios for each compression option, 
and that, if channels are time- varying, current channel state information is available. 
Our algorithm bases decisions purely on this information and does not require a- 
priori knowledge of the packet arrival or channel state probabilities. The algorithm 
is simple to implement and is robust to situations where these probabilities can 
change. This work is important as it demonstrates a principled method of making 
on-line compression decisions in a stochastic system with correlated data. Our 
solution applies the techniques of Lyapunov optimization developed in our previous 
work [Neely 2006] [Georgiadis et al. 2006], and is perhaps the first application of 
these techniques to the dynamic compression problem. This paper also extends the 
general theory by introducing a novel concept of place-holder bits to improve delay 
in stochastic networks with costs. Related Lyapunov optimization techniques for 
network flow control applications are developed in [Neely 2003] [Neely et al. 2005], 
and alternative fluid model approaches are developed in [Stolyar 2005] [Eryilmaz 
and Srikant 2005]. 

In the next section we describe the system model, and in Section 3 we characterize 
the minimum average power in terms of an optimization problem based on channel 
and packet arrival probabilities. In Section 4 we develop an on-line algorithm that 
makes simple decisions based only on current information. The algorithm achieves 
time average power that can be pushed arbitrarily close to optimum via a simple 
control parameter that also affects an average delay tradeoff. A simple improvement 
via place-holder bits is developed in Section 5. Extensions to systems with distortion 
constraints are given in Section 6. Simulations are provided in Section 7. 

2. SYSTEM MODEL 

Consider the wireless link of Fig. 1 that operates in slotted time and receives packets 
from A^ different sensor units. If an individual sensor sends data during a timeslot, 
this data is in the form of a fixed length packet of size b bits, containing sensed 
information. Let A{t) represent the number of sensors that send packets during slot 
t, so that A(t) e {0, 1, . . . , N}. The data from these A{t) packets may be correlated, 
and hence it may be possible to compress the information within the A{t) packets 
(consisting of A{t)b bits) into a smaller data unit for transmission over the wireless 
link. This is done via a compression function ^(a, k) defined as follows. There are 
K +1 compression options, comprising a set /C = {0, 1, . . . , K}. Option represents 
no attempted compression, and options {1, 2, . . . , K} represent various alternative 
methods to compress the data. The function ^(a, k) takes input o G {0,1,..., N} 
(representing the number of newly arriving packets) and compression option /c G /C, 
and generates a random variable output R, representing the total size of the data 
after compression. 

Every timeslot the link controller observes the random number of new packet 
arrivals A{t) and chooses a compression option k{t) G /C, yielding the random com- 
pressed output R(t) — 'i'(A(t),k{t)). Let Pcomp{t) represent the power expended 
by this compression operation, and assume this is also a random function of the 
number of packets compressed and the compression option. We assume that the 
compressed output R{t) is conditionally i.i.d. over all slots that have the same 
number of packet arrivals A{t) and the same compression decision k(t). Likewise, 



compression power Pcomp{t) is conditionally i.i.d. over all slots with the same A{t) 
and k{t). The average compressed output m{a,k) and the average power expendi- 
ture 0(a, k) associated with A{t) — a, k{t) — k are defined: 



m{a, fc) = E {*(A(t), k{t)) \ A{t) = a, k(t) ^ k} 
(j){a,k) = EjPcompW I A{t) = a,k{t) ^ k} 



(1) 
(2) 



We assume the values of m{a, k) and (/)(a, k) are known so that the following table 
can be constructed: 



k 


*{a,fc) 


£{-^(0, fc)} 


E{Pcomp| a,fc} 





ab 


ab 


0{a,O) = 


1 


Random 


m{a, 1) 


0{a,l) 


2 


Random 


m{a, 2) 


0(a,2) 










K 


Random 


m{a, K) 


0(a,^) 



Note that we assume ^(a, 0) = ab and (/)(a, 0) = 0, as the compression option 
k = does not compress any data and also does not expend any power. We 
further assume that m{a, k) < ab for all a e {0,1,..., TV} and all A; G /C, so that 
compression is not expected to expand the data. 

2.1 Data Transmission and Queueing 

The compressed data R{t) = 'i'(A{t),k{t)) is delivered to a queueing buffer for 
transmission over the wireless link (sec Fig. 1). Let U{t) represent the current 
number of bits (or unfinished work) in the queue. The queue backlog evolution is 
given by: 



U{t + 1) = max[t/(t) - i^{t), 0] + R{t) 



(3) 



where fi{t) is the transmission rate offered by the link on slot t. This rate is 
determined by the current channel condition and the current transmission power 
allocation decision, as in [Georgiadis et al. 2006]. Specifically, the channel is as- 
sumed to be constant over the duration of a slot, but can potentially change from 
slot to slot. Let S{t) represent the current channel state, which is assumed to take 
values in some finite set S. We assume the channel state S{t) is known at the 
beginning of each slot t, so that the link can make an opportunistic transmission 
power allocation decision Ptran{t), yielding a transmission rate /i(i) given by: 

^i{t)^C{Ptran{t),S{t)) 

where C'{P, s) is the rate-power curve associated with the modulation and coding 
schemes used for transmission over the channel. We assume C{P, s) is continuous 
in power P for each channel state s G 5. Transmission power allocations P{t) 
are restricted to some compact set V for all slots t, where V contains a maximum 
transmission power Pmax- For example, the set V can contain a discrete set of 
power levels, such as the two element set V — {0, Pmax]- Alternatively, V can be a 
continuous interval, such as "P = {P | < P < Pmax}- We assume throughout that 
e T' and that C(0, s) = for all channel states s S 5, so that zero transmission 
power yields a zero transmission rate. Further, we assume that C{PmaxTs) > 



C{P, s) for all s G iS and all P £ P, so that allocating maximum power yields the 
largest transmission rate that is possible under the given channel state. 

2.2 Stochastic Assumptions and the Control Objective 

For simplicity, we assume the packet arrival process A{t) is i.i.d. over slots with a 
general probability distribution pa{<i) = Pr[A{t) — a]. Likewise, the channel state 
process S{t) is i.i.d. over slots with a general distribution tTs ~ Pr[S{t) = s]} 
The distributions pa (a) and tTs are not necessarily known to the link controller. 
Every slot the link controller observes the number of new packets A{t), the current 
queue backlog U{t), and the current channel state S{t), and makes a compression 
decision fc(t) € /C (expending power Pcompit)) and a transmission power allocation 
Ptran{t) G P- Thc total time average power expenditure is given by: 

1 *"^ 

lim -y^JPcornpiT) + Ptranir)] 
r=0 

The goal is to make compression and transmission decisions to minimize time aver- 
age power while ensuring the queue U{t) is stable. Formally, we define a queucing 
process U{t) to be stable if: 

t-i 



limsup-^E{[/(r)} < 



T=0 

This type of stability is often referred to as strong stability, as it implies a finite 
average backlog and hence a finite average delay. In Section 4, we shall design a 
class of dynamic algorithms that can drive time average power arbitrarily close to 
the minimum average power required for stability, with a corresponding explicit 
tradeoff in average queue backlog and average delay. 
Define rmin and rmax as follows: 

r„„ A E|minm(A(t),fc)| (4) 

r^ax = E{C(P„a,,5(t))} (5) 

where the expectations are taken over the randomness of A{t) and S{t) via the 
distributions Pa{<i) and tt^. Thus, Tmin is the minimum average bit rate delivered 
to the queueing system (in units of bits/slot), assuming the compression option 
that results in the largest expected bit reduction is used every slot. The value Vmax 
represents the maximum possible average transmission rate over the wireless link. 
We assume throughout that rmin < fmax, so that it is possible to stabilize the 
system. 

Thus, there are two reasons to compress data: (i) In order to stabilize the queue, 
we may need to compress (particularly if E {A(t)} b > r^nax)- (ii) We may actually 



-"^ Using thc T-slot Lyapunov drift techniques described in [Georgiadis et ah 2006], our analysis 
can be generaUzed to show that the same algorithms we derive under the i.i.d. assumption yield 
similar performance for arbitrary ergodic arrival and channel processes A(t) and S{i), with delay 
bounds that increase by a constant factor related to thc mixing times of the processes. 



save power if the power used to compress is less than the extra amount of power 
that would be used transmitting the extra data if it were not compressed. 

2.3 Discussion of the System Model 

This simple model captures a wide class of systems where data compression is 
important. The N sensor scenario of Fig. 1 captures the possibility of randomly 
arriving data that is spatially correlated. An example is when there are multiple 
sensors in an environment and only a random subset of them detect a particular 
event. The data provided by these sensors is thus correlated but not necessarily 
identical, as each observation can offer new information. 

The case of compression because of time correlated data can also be treated in this 
model by re-defining N to represent the time over which a frame of data samples 
are gathered. Indeed, suppose a timeslot t is composed of N mini-slots, where 
data can arrive on any or all of the mini-slots. The value of A{t) now represents 
the random number of packets arriving over the N mini-slots, and the compression 
functions m(a, k) and (/)(a, k) now represent averages associated with compressing 
the timc-corrclatcd data. This of course assumes compression is contained to data 
arriving within the same frame, and does not treat inter-frame compression. 

Our time- varying channel model is useful for systems with mobility, environmen- 
tal changes, or restrictions that create time-varying transmission opportunities. 
This allows for opportunistic scheduling which can help to further reduce power 
expenditure. We do not consider the additional power required to measure the 
channel conditions here. Extensions that treat this issue can likely be obtained 
using the techniques for optimizing measurement decisions developed in [Li and 
Neely 2007]. A special case of the time- varying channel model is the static channel 
assumption, where S{t) is the same for all timeslots t. This special case is similar 
to the static assumption in [Barr and Asanovic 2003] [Sadler and Martonosi 2006] . 
However, this static channel scenario still creates an interesting problem that is 
much different from [Barr and Asanovic 2003] [Sadler and Martonosi 2006]. In- 
deed, the random packet arrivals (with raw data rate that is possibly larger than 
link capacity) and the potentially non-linear rate-power curve necessitate a dynamic 
compression strategy that is not obvious, that depends on the packet arrival dis- 
tribution, and that docs not necessarily use the same compression option on every 
slot. 

Here we assume that the compression options available within the set K, are 
sufficient to ensure that the resulting data transmitted over the link has an ac- 
ceptable fidelity. An example is lossless data compression, such as Huffman or 
Lempel-Ziv source coding, where all original data packets can be reconstructed at 
the destination. Alternatively, we might have some compression options k G K, rep- 
resenting lossy compression, provided that the distortion that may be introduced 
is acceptable. Extensions to systems that explicitly consider distortion due to lossy 
compression are considered in Section 6. 

3. MINIMUM AVERAGE POWER 

Here we characterize the minimum time average power required for queue stability. 
We first define separate functions h*{r) and g*{r) that describe the minimum av- 
erage power for compression and transmission, respectively, over a restricted class 



of stationary randomized algorithms. These functions depend on the steady state 
arrival and channel distributions pA (a) and tt^ . We then show that these functions 
can be used to define system optimality over the class of all possible decision strate- 
gies, including strategies that do not necessarily make stationary and randomized 
decisions. 

3.1 The Functions h*{r) and g*{r) 

Definition 1. For any value r such that rmin < ^ < UE,{A{t)}, the minimum- 
power compression function h* (r) is defined as the infimum value h for which there 
exist probabilities {'^a,k) for a 6 {0, 1, . . . , N}, k € IC, such that the following con- 
straints are satisfied: 

N K 

'^^PA{a)la,k(t){a, k) ^ h (6) 

a=Ofe=l 
N K 

^^PA{a)la.km{a, k) < r (7) 

a=Ofe=l 

7a. fe > for all a, k (8) 

K 
y^ 7a. fe = 1 for all a (9) 

fe=i 

Intuitively, the (7a, fe) values define a stationary randomized policy that observes 
the current arrivals A{t) and uses compression option k with probability 7a, fe when- 
ever A{t) = a. The expression on the left hand side of (6) is the expected compres- 
sion power E {Pcomp{t)} for this policy. Likewise, the expression on the left hand 
side of (7) is the expected number of bits E {R{t)} at the output of the compressor 
for this policy. The value of h* (r) is thus the smallest possible average power due 
to compression, infimized over all such stationary randomized policies that yield 
E{i?(t)} < r. Note from (4) that it is possible to have a stationary randomized 
policy that yields E {R{t)} = Vmim and hence the function h*{r) is well defined for 
any r > rmin- Further, the following lemma shows that the infimum value h*{r) 
can be achieved by a particular stationary randomized algorithm. 

Lemma 1. For any r such that rmin < ?' < bE{A(t)}, there exists a particular 
stationary randomized policy that makes compression decisions k*{t) as a random 
function of the observed A{t) value (and independent of queue backlog), such that: 

E{0(A(i),r(i))} = h*{r) (10) 

E{m{A{t),k*{t))} ^ r (11) 

where the above expectations are taken with respect to the steady state packet arrival 
distribution pa (a) and the randomized compression decisions k* (t) . 

Proof. The proof follows by continuity of the functions on the left hand side of 
(6) and (7) with respect to 7a, fe, and by compactness of the set of all (7a, fe) that 
satisfy (8) and (9). See Appendix A for details. D 

Similar to the function h*(r), we define g*{r) as the smallest possible average 
transmission power required for a stationary randomized algorithm to support a 



transmission rate of at least r. The precise definition is given below. 

Definition 2. For any value r such that < r < Tmax, the minimum-power 
transmission function g* (r) is defined as the infimum value g for which there exists 
a stationary randomized power allocation policy that chooses transmission power 
Ptran(t) as a random function of the observed channel state S{t) (and independent 
of current queue backlog), such that: 

E{Ptran{t)} - 5 (12) 

E{C(Ptran{t),S{t))} > r (13) 

The function g*{r) is well defined whenever r < r^ax because it is possible to 
satisfy the constraint (13). Indeed, note by (5) that the policy Ptran{t) = Pmax 
for all t yields E{C(Ptran(i), 'S'(i))} = Tmax- Furthermore, it is easy to show that 
the inequality constraint in (13) can be replaced by an equality constraint, as any 
policy with an average transmission rate larger than r can be modified to achieve 
rate r exactly while using strictly less power. This can be done by independently 
setting Ptran{t) ~ with some probability every slot, yielding a zero transmission 
rate in that slot. 

Because the set V is compact and the function C{P, s) is continuous in power P 
for all channel states s e 5, an argument similar to the proof of Lemma 1 can be 
used to show that the infimum average power g*{r) can be achieved by a particular 
stationary randomized policy^ Specifically, for any r such that < r < r^ax, there 
exists a stationary randomized algorithm that chooses transmission power Pf^an (t) 
that yields: 

E{C(P,;„„(i),5(t))} = r (14) 

^{Kanit)} - 9*{r) (15) 

3.2 Structural Properties of /i*(r) and g*(r) 

It is not difficult to show that h*(r) is a non- increasing function of r (because 
less compression power is required if a larger compressor output rate is allowed), 
and that g*{r) is a non-decreasing function of r (because more transmission power 
is required to support a larger transmission rate). Further, both functions are 
convex. It is interesting to note that in the special case when there is no channel 
state variation so that C{P,s) = C{P), and when the function C{P) is strictly 
increasing and concave, then g*{r) = C~^{r), i.e., it is given by the inverse of 
C{P). Details on the structure of the g*(r) function in the general time- varying 
case are given in [Neely 2006]. 

3.3 Minimum Average Power for Stability 

The following theorem establishes the minimum time average power required for 
queue stability in terms of the h*(r) and g*{r) functions. We consider all possible 
algorithms for making compression decisions k(t) G /C and transmission power 



^More generally, the infimum can be achieved whenever C{P, s) is upper semi-continuous in P for 
every channel state s £ <S (see [Bertsekas et al. 2003] for a definition). The upper semi-continuity 
property is a mild property that is true of every practical curve C{P, s). All results of this paper 
hold when continuity is replaced by upper semi-continuity. 



decisions Ptran{t) G V over time, including algorithms that are not necessarily in 
the class of stationary randomized policies. 

Theorem 1. Let A{t) and S{t) he ergodic with steady state distributions pA^a) 
and -Kg, respectively (such as processes that are i.i.d. over slots, or more general 
Markov modulated processes). Assume that rmin < f max (defined in (4), (5)). Then 
any joint compression and transmission rate scheduling algorithm that stabilizes the 
queue U{t) yields a time average power expenditure that satisfies: 

t-i 



limSUp - ^ E {Pcornpir) + Ptranir)} > P* 
t — >oo t 



T=0 

where P*^ is defined as the optimal solution to the following problem: 

Minimize: h*{r)+g*{r) (16) 

Subject to: r^in <r< mm[rrnax , blK {A{t)}] (17) 

Proof. See Appendix B. D 

The above theorem shows that time average power must be greater than or equal 
to P*y for queue stability. The result can be understood intuitively by observing 
that if r is the rate of bits arriving to the queue from the compressor, then average 
transmission power can be minimized while maintaining stability by pushing the 
time average transmission rate down closer and closer to r. The optimization 
problem corresponding to this definition of P*^ may be difficult to solve in practice, 
as it would require exact knowledge of the h*{r) and g*{r) functions, which in 
turn requires full a-priori knowledge of the distributions pa (a) and tt^ . In the next 
section, we design a simple class of dynamic algorithms that stabilize the queue 
without this knowledge, and that push time average power arbitrarily close to P*^. 

4. THE DYNAMIC COMPRESSION ALGORITHM 

Our dynamic algorithm is decoupled into separate policies for data compression and 
transmission rate scheduling. It is defined in terms of a control parameter V > 
that affects an energy-delay tradeoff. 

The Dynamic Compression and Transmission Algorithm: 

Compression: Every slot t, observe the number of new packet arrivals A(t) and 
the current queue backlog U{t), and choose compression option k{t) £ /C as follows: 

k{t) = argmin [U {t)m{A{t) , k) + F0(A(t), k)] (18) 

If there are multiple compression options fc G /C that minimize U (t)m{A{t) , k) + 
V(j){A{t), k), break ties arbitrarily. 

Transmission : Every slot t, observe the current channel state S{t) and the current 
queue backlog C/(i), and choose transmission power Ptran{t) 6 "P as follows: 

Ptran{t) = arg max [U{t)C{P, S{t)) - VP] (19) 

Recall that V is assumed to be compact and the C{P^ s) function is upper semi- 
continuous, and hence there exists a maximizing power allocation. If there are 
multiple power options that maximize U(t)C{P, S{t)) — VP, break ties arbitrarily. 
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The compression policy involves a simple comparison oi K + 1 values found 
by evaluating the m(a, k) and 0(a, k) functions for all k E IC, and can easily be 
accomplished in real time. The transmission policy is a special case of the Energy 
Efficient Control Algorithm (EECA) policy developed in [Neely 2006] , and typically 
can also be solved very simply in real time. The next theorem establishes the 
performance of the combined algorithm. 

Theorem 2. (Algorithm Performance) Suppose packet arrivals A{t) are i.i.d. 
over slots with distribution pA{a), and channel states S{t) and are i.i.d. over slots 
with distribution tt^. For any control parameter V > 0, the dynamic compression 
and transmission scheduling algorithm yields power expenditure and queue backlog 
that satisfy the following: 

Plot < PL + B/V (20) 

where Ptot o,nd U are the time averages for power expenditure and queue backlog, 
defined: 

- 1 *"^ 

Ptot = limSUp - V" E {Pcomp{T) + Ptran{T)} 



*-- ^ r=0 



t ^—' 

T = 



V 4 limsup-2^E{C/(T)} 
and where B and 4'max o,fe constants given by: 



B A ^[a^+E {C{P^ax,S{t)f}] (22) 



max[</)(A(t),fc)]| (23) 

where a^ is an upper bound on E |i?(t)^} for all slots t. For example, if no com- 
pression operation expands the data, then R{t) — '^{A{t),k) < bA{t) for all t, and 
hence a^ is defined: 

a2Afe2E{A(t)2} 

We prove Theorem 2 in the next subsection. Note that the parameter y > 
can be chosen to make B jY arbitrarily small, ensuring by (20) that time average 
power is arbitrarily close to the optimal value P*^. However, the resulting average 
queue backlog bound grows linearly with Y . By Little's Theorem, the average 
queue backlog is proportional to average delay [Bertsekas and Gallager 1992]. This 
establishes an explicit tradeoff between average power expenditure and delay. 

As an implementation detail, we note for simplicity that we can use units of 
bits, bits/slot, and milli- Watts for U(t), C{P,S), and P. However, these units are 
arbitrary and any consistent units will work, with performance given by (20) and 
(21). Indeed, any unit changes are captured in the Y constant (where Y has units of 
bits /mW for the units above). For example, if milli- Watts are changed to Watts, 
then the algorithm will make the exact same control decisions for k{t) and Ptran{t) 
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over time, and hence yields the exact same sample path of energy use and queue 
backlog, as long as the V constant is appropriately changed by a factor of 1000. If 
bits are changed to kilobits, then V must change by a factor of 10^. 

4.1 Lyapunov Performance Analysis for Theorem 2 

Our proof of Theorem 2 relies on the performance optimal Lyapunov scheduling 
techniques from [Georgiadis et al. 2006] [Neely 2006]. First define the following 
quadratic Lyapunov function of queue backlog U{t): 

L(C/(i))Ai[/(i)2 

Define the one-step conditional Lyapunov drift A(U{t)) as follows:^ 

A(L/(t)) Ae {LiUit + 1)) - LiUit)) I Um (24) 

The following simple lemma from [Georgiadis et al. 2006] shall be useful. 

Lemma 2. (Lyapunov drift [Georgiadis et al. 2006]) Let L{U{t)) be a non-negative 
function ofU{t) with Lyapunov drift A(C/(i)) defined in (24)- If there are stochastic 
processes a{t) and f3{t) such that every slot t and for all possible values ofU{t), the 
conditional Lyapunov drift satisfies: 

A(C/(t)) < E {I3{t) - a{t) I U{t)} (25) 

then: 

t-i -. t-i 

limsup-^E{a(T)} < limsup - ^ E {/3(t)} 

The proof involves taking expectations of (25), using iterated expectations, and 
summing the resulting telescoping series (see [Georgiadis et al. 2006] for details). 

The queue backlog U{t) for our system satisfies the queue evolution equation (3). 
Specifically, the queue has arrival process R{t) = 'i'(A{t),k{t)) and transmission 
rate process /i(i) — C{Ptran{t), S{t)), where the k{t) and Ptran{t) control decisions 
are determined by the dynamic compression and transmission algorithm of the 
previous sub-section. The Lyapunov drift is given by the following lemma. 

Lemma 3. (Computing A{U{t))) Under the queue evolution equation (3) and us- 
ing the quadratic Lyapunov function L{U{t)) = ■^U{t)'^ , the Lyapunov drift A{U{t)) 
satisfies the following for all t and all U{t): 

A{U{t)) <B- U{i)¥. {fi{t) - m{A{t), k{t)) \ U{t)} (26) 

where /i(t) = C{Ptran{t), S{t)), and B is given in (22). The expectation above 
is taken with respect to the random channels and arrivals S(t) and A(t), and the 
potentially random control actions k{t) and Ptran{t). 



^More complete notation would be A(f7(t),t), as the drift depends on the scheduling policy 
which may also depend on time t. However, wc use the simpler notation A(f/(t)) as a formal 
representation of the right hand side of (24). See [Georgiadis et al. 2006] for further details on 
Lyapunov drift. 
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Proof. From (3) we have: 

lu(t + lf = i(max[t/(t)-M(t),0]+i?(t))' 

< ^[Uitf+fi{tf + Ritf] 
~U{t){^i{t) - R{t)) 
and hence (taking conditional expectations given U{t)): 

AiUit)) < ^E{,4tf + R{tr\Uit)} 

-U{t)E{^{t)-R{t)\ U{t)} 

It is clear that the value ^E |/i(t)^ + i?(t)^ | U{t)^ is less than or equal to the 
constant B defined in (22), and hence: 

A{U{t)) <B- U{t)E {^{t) - R{t) I U{t)} (27) 

Noting that R{t) = ^(A(t), k{t)) and using iterated expectations, we have: 

E{R{t)\ U{t)} 

= E{^{A{t),k{t))\U{t)} 

= E{E{^{A{t),k{t)) I U{t),A{t),k{t)} I U{t)} 

= E{m{A{t),k{t))\ U{t)} 

where we have used the definition of m{a,k) given in (1). Using this equality in 
(27) yields the resuh. D 

Following the Lyapunov optimization framework of [Georgiadis et al. 2006] [Neely 
2006], we add a weighted cost term to the drift expression. Specifically, from (26) 
we have: 

A(C/(t)) + VE{Pcomp{t) + Ptran{t) \ U {t)) < 

B - U(t)E {C{Ptran{t). S{t)) - m{A{t),k{t)) I U{t)} 

+VE {Pcompit) + Ptranit) \ U (t)} (28) 

where we have just added an additional term to both sides of (26). Note that 
^ {Pcomp{t) I U{t)} can be expressed as follows (using iterated expectations): 

E{Pcomp(i)| U{t)} 

= E{E{P,omp{t) 1 U{t),A{t),k{t)} \ U{t)} 
= E{cb{A{t),k{t))\U{t)} 

Using this equality in the right hand side of (28) and re-arranging terms yields: 

A{U{t)) + VE{P,orr.p{t) + Ptranit) \ U (t)} < 

B-E {U{t)C (Ptranit), S{t)) - V Ptran{t) \ U [t)] 

+E{U{t)m{A{t),k{t)) + V(t){A{t),k{t))\ U{t)} (29) 

Now note that we have not yet used the properties of the dynamic compression 
and transmission policy. Indeed, the above expression (29) is a bound that holds 
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for any compression and transmission scheduling decisions k{t) G /C, Ptran{t) G "P 
that are made on slot i, including randomized decisions. However, note that the 
dynamic compression and transmission strategy is designed specifically to minimize 
the right hand side of (29) over all alternative decisions that can be made on slot 
t. Indeed, the compression algorithm observes A{t) and U{t) and chooses k{t) e JC 
to minimize U{t)m{A{t), k{t)) + V(j){A{t), k{t)), which thus minimizes the following 
term over all alternative decisions that can be made on slot t: 

E {U{t)m{A{t),k{t)) + V(j){A{t), k{t)) I U{t)} 

Similarly, the transmission power allocation algorithm is designed to minimize the 
following term over all alternative decisions that can be made on slot t: 

-E{U{t)C{Ptran{t),S{t)) - V Ptran{t) \ U (t)} 

It follows that the right hand side of (29) is less than or equal to the corresponding 
expression with Ptran{t) and k{t) replaced by Pfrani^) and fc*(i), where Pt^ani^) 
and k*{t) are any other (possibly randomized) policies that satisfy Ptrani^) ^ ^ 
and k*{t) G /C: 

A{U{t)) + VE{P,omp{t) + Ptranit) \ U {t)} < 

B -E{u{t)c{p:,,„{t), s{t)) - yp;™„(t) I u{t)} 

+E{U{t)m{A{t),k*{t)) + V<l){A{t),k*{t))\ U{t)} (30) 

Now let ri be any particular value that satisfies rmin < ?'i < bK{A{t)}, and let 
k*{t) be the stationary randomized policy that yields: 

E{0(A(t),fc*(i))} = h*{n) (31) 

E{m(A(i),fc*(t))} = ri (32) 

Such a policy exists by (10) and (11) of Lemma 1. Similarly, let r2 be any value 
that satisfies < r2 < r^nax, and let Pt^ani^) ^'e the stationary randomized power 
allocation policy that yields: 

¥.{C{P;,,^{t),S{t))} = r2 (33) 

^ {Ptranit)} = 5*('^2) (34) 

Such a policy exists by (14) and (15). Further, the stationary randomized policies 
of (31)-(34) base decisions only on the current A{t) and S{t) states, which are i.i.d. 
over slots (and are hence independent of the current queue backlog U(t)). Thus, 
the expectations of (31)-(34) are the same when conditioned on U{t). Plugging 
(31)-(34) into the right hand side of (30) thus yields: 

A(C/(i)) + VE{Pcornp{t) + Ptranit) \ U [t)} < 

B - U{t){r2 - ri) + V{h*{n) + g^r^)) (35) 

The above inequality holds for all ri and r2 that satisfy rmin < ^i < 6E{A(t)} 
and < r2 < Tmax- Let ri = r2 = r* , where r* is the value of r that optimizes the 
problem in (16) and (17) of Theorem 1, so that P*^ = h*(r*) +g*{r*). Plugging 
into (35), we have: 

HU{t)) + FE {P,omp{t) + Ptranit) \ U {t)} <B + V P*^^ 
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Using the above drift inequality in the Lyapunov Drift Lemma (Lemma 2) and 
defining a{t) = VPcompit) + VPtran{t) and /3(t) = S + FP,*„ yields Ptot < PL + 
B/V, proving equation (20) of Theorem 2. 

Now choose ri = rmin and r2 = Vmax- Plugging into (35) and noting that 
Pcompit) > and Ptran{t) > gives: 

A{Uit)) <B- U{t){rraax ~ r„„„) 
+ V{h*{r„iin) +g*{rrnax)) 

Using the above drift inequality in the Lyapunov Drift Lemma (Lemma 2) and 
defining a{t) = U(t){rmax ~ rmin) and /3(t) = B + V{h*{rmin) + g*{rmax)), we 
have: 

*"' B + Vih*ir„nn)+9*{rmax)) 



limsup-^E{C/(T)} < 



t — *oo I- _j-, [rmax ^jYiin) 



The result of (21) follows because g*{rmax) < Pmax and h*{rmin) < (pmax- This 
completes the proof of Theorem 2. 

5. A SIMPLE DELAY IMPROVEMENT 

Here we present a simple improvement to the transmission algorithm that can de- 
crease queue backlog while maintaining the exact same average power performance 
specified in Theorem 2. First observe that the performance theorem (Theorem 2) 
and the Lyapunov Drift Lemma (Lemma 2) both specify time average behavior that 
is independent of the initial queue backlog. Indeed, the affects of the initial con- 
dition are transient and decay over time. Now suppose the transmission algorithm 
has the following property: 

Property 1: There exists a finite constant Uthresh > such that if t/(0) > Uthresh, 
then U(t) > Uthresh for all time t > 0. 

Thus, Property 1 says that if the queue has an initial condition of at least Uthresh 
bits, then it will never fall below this threshold of bits. Clearly Property 1 always 
holds with Uthresh = 0. In any system where Property 1 holds for some constant 
Uthresh > 0, thcu thc first Uthresh bits in the queue are just acting as a place holder 
to make U{t) large enough to properly affect the stochastic optimization. 

5.1 Example Showing that Uthresh > is Typical 

Suppose there is a finite constant (imax such that: 

C(P, 5*) < PrnaxP for idl P £ V and all S* e 5 

For example, if the transmission rate function C(P, 5*) is differentiable with respect 
to P, then /3max can be defined as the largest derivative with respect to P over all 
possible channel states. Because the algorithm chooses Ptran(t) every slot as the 
maximizer of C/(i)C(P, £'(<)) - VP over all P e P, it is clear that Ptran{t) = 
whenever U{t)Pmax < V. Indeed, we have for any channel state S{t): 

U{t)C{P, S{t)) -VP< [U{t)Pmax - V]P 

which is maximized only by P = if U(t)f3max < V. Therefore, if U{t) < V/ Pmax, 
we have Ptran{t) = and hence ii{t) = C(Ptran{t), S(t)) = 0, so that thc queue 
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backlog cannot further decrease. Define Umax as the largest possible transmission 
rate during a single slot (equal to the maximum of C{Pmax, S) over all S* G 5). It 
follows that Property 1 holds in this example with:^ 



Uthresh^T^a^ 



V 
0, -5 M>i 

Pmax 



(36) 

The value Uthresh determines the number of place holder bits required in the system. 

5.2 Delay Improvement Via Place Holder Bits 

If Property 1 holds for Uthresh > 0, performance can be improved in the following 
way: With U(t) being the actual queue backlog, define the place-holder backlog 
U{t) as follows: 

U{t)^U{t) + Uthresh 

The value U{t) can be viewed as backlog that is equal to the actual backlog plus 
Uthresh "fake bits." Now assume that U{0) — 0, but implement the Dynamic 
Compression and Transmission Algorithm using the place-holder backlog U{t) ev- 
erywhere, instead of the actual queue backlog. That is, choose k{t) £ K, and 
Ptran{t) G P as foUows: 

k{t) = aigminp {t)m{A{t),k) + V(t>{A{t),k)] 

Ptranit) - Sirgm^x[1J {t)C{P, S (t)) ^ V P] 

With this implementation, we have U{0) = Uthresh, and so U{t) > Uthresh for all 
t (by Property 1). It follows that any transmission decisions never take U{t) lower 
than Uthresh, which is equivalent to saying that all transmission decisions transmit 
only actual data (so that U{t) > 0), rather than fake data. The resulting deci- 
sions are the same as those in a system with initial backlog of Uthresh, which yields 
the same 0{1/V) energy performance guarantee as before, and yields the same 
0{V) time average backlog guarantee for the U{t) backlog. However, the actual 
queue backlog is exactly Uthresh bits lower than tJ{t) at every instant of time, and 
so the time average backlog is also exactly Uthresh bits lower. Thus, this simple 
improvement yields less actual queue backlog in the system, without any loss in 
performance. This improvement does not change the [0(1/F), 0(\^)] tradeoff rela- 
tion (it simply multiphes the 0{V) congestion bound by a smaller coefficient), but 
can yield practical backlog and delay improvements for implementation purposes. 

6. DISTORTION CONSTRAINED DATA COMPRESSION 

In the previous sections, we have assumed that all compression options fc G /C are 
either lossless, or that the distortion they introduce is acceptable. In this section, 
we expand the model to allow each compression option to have its own distortion 
properties. Specifically, let D{t) be a non-negative real number that represents a 
measure of the distortion introduced at time t due to compression. We assume 



^The Uthresh value in (36) satisfies Property 1, but is not necessarily the largest value that 
satisfies this property. Performance can be improved if a larger value that satisfies Property 1 can 
be found, which is often possible for concave ratc-powrcr curves defined over a continuous interval. 
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that this is a random function of A{t) (the number of packets compressed) and 
k{t) (the compression option), and that this function is stationary and independent 
over all slots with the same A{t) and k{t). Define d{a, k) as the expected distortion 
function: 

d{a, k)^E {D{t) I A{t) = a, k{t) = k} 

We assume the maximum second moment of distortion is bounded by some finite 
constant 5"^: 

S^A max [E {D(t)^ I A{t) = a, kit) = k]] 

ae{l,...M}.kelC L I- ^ ' I ^ ■' JJ 

Assuming that distortion is additive, the goal in this section is to make joint 
transmission and compression actions to minimize time average power expenditure 
subject to queue stability and subject to the constraint that time average distortion 
is bounded by a constant day'- 

1 *"^ 
\im -Ye {D{T)}<dav (37) 

T = 

Let )C now represent a set of extended compression options, which includes lossy 
compression options (and may also include the maximum distortion option of throw- 
ing away all data that arrives on slot t). Lossless compression options may still be 
available and yield D{t) = 0. 

To ensure the distortion constraint (37) is satisfied, we introduce a distortion 
queue X{t) that accumulates the total amount of distortion in excess of dav This 
is similar to the virtual power queue introduced in [Neely 2006] for ensuring aver- 
age power constraints, and is an example of a virtual cost queue from [Georgiadis 
et al. 2006]. Specifically, the X{t) queue is implemented purely in software. It is 
initialized to X{0) = 0, and is changed from slot to slot according to the following 
dynamics: 

X{t+l) = iaaii[X{t)-dav,0]+D{t) (38) 

where D{t) is the random amount of distortion introduced on slot t (observed at 
the end of the compression operation). Stabilizing the X{t) queue ensures the time 
average rate of the D{t) "arrivals" is less than or equal to dav, which is equivalent 
to the distortion constraint (37). 

6.1 The Distortion-Constrained Algorithm 

The queue backlog U{t) evolves as before, according to the dynamics (3). The 
transmission algorithm that selects Ptran{t) £ P every slot is the same as before 
(equation (19)), and the new compression algorithm is given as follows: 

Distortion- Constrained Compression Algorithm: Every slot t, observe the num- 
ber of new packet arrivals A(t) and the current queue backlogs U{t) and X{t), and 
choose compression option fc(t) e /C as follows: 

k{t) ^ argmin ^{t)m{A{t),k) + X{t)d{A{t), k) 

kGtC 

+F0(A(t),fc)] 
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After the compression operation, observe the actual distortion D{t) and update 
X{t) according to (38). 

We note that the distortion-constrained dynamic compression algorithm above is 
useful even in cases when the average power expenditure 4'{A{t), k) due to compres- 
sion is negligible. This is because compression can significantly save transmission 
power, although intelligent compression strategies are required to meet the distor- 
tion constraints. 

6.2 Minimum Average Power with Distortion Constraints 

Define h'^ir) as the infimum time average power over all compression strategies 
that make decisions k(t) € IC as a stationary and random function of the observed 
number of packets A{t), subject to a time average output rate of the compressor 
that is at most r bits/slot, and subject to a time average distortion of at most 
dav Let rd,min represent the smallest possible time average output rate of the 
compressor that yields a time average distortion of at most dav , optimized over all 
algorithms that choose k{t) £ /C as a stationary and random function of A{t) (and 
hence independently of queue backlog). We assume that rd,min < fmax, so that it 
is feasible to meet the distortion constraint. We further assume that it is possible 
to meet the distortion constraint with strict inequality while stabilizing the system. 
That is, we assume there is an e > such that it is possible to achieve a time 
average distortion rate oi D < dav ^ e using a policy that chooses k{t) G /C as a 
stationary and random function of A{t), such that E {m{A{t)^ ^(0)} ^ fmax- 

Recall that g* (r) is the infimum time average transmission power over all station- 
ary randomized strategies Ptran(t) G V that yield time average transmission rate 
at least r bits/slot. Let P^„ represent the minimum time average power (summed 
over compression and transmission powers) required to stabilize the queueing sys- 
tem subject to the distortion constraint (37). 

Theorem 3. (Distortion Constrained Minimum Average Power) The distortion- 
constrained minimum time average power P*^ is equal to the solution of the follow- 
ing optimization problem: 

Minimize: h'^(r)+g*{r) (39) 

Subject to: rd,min <r< Tam[rmax,bE {A{t)}] 

Proof. The proof is similar to the proof of Theorem 1 and is omitted for 
brevity. D 

6.3 Lyapunov Analysis 

Let Z{t)=[U{t), X{t)] be the combined queue state, and define the Lyapunov func- 
tion: 

L(Z(t))Alt/(t)2 + Ixitf 
Define the Lyapunov drift: 

A(Z(t)) Ae {LiZit + 1)) - LiZit)) I Z(t)} 
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Lemma 4. For any constant V > 0, the Lyapunov drift A{Z{t)) satisfies the 
following for all t and all Z{t): 

/^{Z{t)) + V¥.{Ptot{t)\Z{t)}<C 

-U{t)¥.{C{Ptran{t).S{t))\ Z{t)} 

+U{t)¥.{m{A{t),k{t))\ Z{t)] 

-X(t)E {dar, - d{A{t),k{t)) I Zit)} 
+ VE {Ptranit) + 0(A(i), Ht)) \ Z (t)} 

where the constant C is given by: 

C ^l K. + -5' + ^2 + E{C{Ema.,S{t))^}] (40) 

Proof. The proof is similar to the proof of Lemma 3 and is omitted for brevity. D 

It can be seen that the Distortion Constrained Compression and Transmission 
Algorithm is designed to observe current queue backlogs X{t), U{t) and arrival 
and channel states A{t), S{t), and take control actions k{t) G /C, Ptran(t) G "P to 
minimize the right hand side of the drift bound given in Lemma 4. 

Theorem 4. (Algorithm Performance with Distortion Constraint) Suppose A{t) 
and S{t) are i.i.d. over slots, and that rd,min < fmax- For any control parame- 
ter V > Q, the Distortion Constrained Compression and Transmission Algorithm 
stabilizes the network and satisfies: 

Ptot < PL + f (41) 

Tj ^ ^ K^rnax + (Pmax) ( A0\ 

[J'max '^d,min) 
D < dav (43) 

where U,Ptot, o,nd D represent limsup time average expected queue backlog, total 
power, and distortion, and the constant C is defined in (40). 

Proof. See Appendix D. D 

Thus, the algorithm meets the time average distortion constraint, and the param- 
eter V can be used to push total time average power expenditure within 0{1/V) of 
the optimal P*,^, with an 0{V) tradeoff in average queue congestion U and hence 
average delay. We note that an improved delay performance can be achieved by 
using U{t) = U{t) + Uthresh as a replacement for U(t), as described in Section 5, 
with Uthresh Satisfying Property 1, such as the value given in (36). 

7. SIMULATIONS 

For simplicity, we consider simulations of the dynamic compression and transmis- 
sion algorithm of Section 4 (with the simple improvement of Section 5), without 
treating distortion constraints. To begin, we first consider a system where the 
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optimal compression decision is trivial and does not require a stochastic optimiza- 
tion. Specifically, suppose that we have a system where all three of the following 
"Singularity Assumptions" hold: 

— The channel is static, so that S{t) is the same for all t. 

— The rate-power curve is linear in power, so that C{P) = aP for all P G P, for 
some constant a. 

— The raw data arrival rate is less than the maximum transmission rate, that is, 
VK{A{t)}<aP^ax- 

In this simple case, the time average transmission power is directly proportional to 
the time average rate of bits transmitted, and so we do not require careful transmis- 
sion decisions (all transmissions are equally energy efRcient). Further, compression 
is not required for stability. It is easy to show in this special case that the exact 
minimum energy expenditure is achieved by the alternative algorithm of observing 
A{t) every slot t and choosing a compression option k{t) e /C as follows: 



k{t) = argmin 



(p(A(t),k) + -m(A(t),k) 
a 



(44) 



and then transmitting whenever there is a sufficient amount of backlog to achieve 
an efhciency of a bits/unit energy. That is, we simply choose the compression op- 
tion that minimizes the sum of the total energy required to compress and transmit 
the bits. A similar observation is used in [Barr and Asanovic 2003] in the study 
of compression energy ratios for popular algorithms. However, this k{t) policy is 
fragile, in that its optimality strongly relies on all three of the above Singularity As- 
sumptions. Our dynamic compression and transmission algorithm is an all-purpose 
algorithm that should work well for any system, including systems that satisfy the 
above Singularity Assumptions, as well as systems that do not. 

7.1 Scenario I: Singularity Assumptions 

We first consider a scenario where all three of the "Singularity Assumptions" hold. 
The channel is static with S{t) = ON for all time slots t. The transmit power is 
constrained to two options V = {0,1} (we used normalized units of power). The 
rate-power curve is given by C{P = 1) = 2048 bits/slot, and C{P = 0) = 0. It 
is clear that the optimal transmission decision in this scenario is to transmit only 
when the queue size is greater than or equal to fimax — 2048 bits, so that all 
transmissions have efficiency a — fimax bits per unit power. 

The wireless link receives data from 8 different sensor units. The packet arrival 
process at each sensor is i.i.d. over slots and follows a Bernoulli distribution with 
the probability of an arrival p = ^. We fix the packet size, b, to 256 bits. Hence, 
the arrival process for the wireless link, A{t), follows a Binomial distribution with 
parameters (8, ^) with an average arrival rate of ML{A{t)} = 1024 bits/slot. Note 
that in this case we have bK{A{t)} < Umax, and so compression is not needed for 
queue stability. 

Two compression options are available to the link controller (/C = {0,1}). For 
A{t) = a and k{t) = 1, the size of the data after compression, R{t), is uniformly 
distributed in [^,%^] and the compression power is uniformly distributed in 
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Fig. 2. Avg. Power expenditure vs. V (two nearly identical curves are shown). 
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Fig. 3. Avg. Queue size vs. V 



[0.45,0.55]. Hence, the average compressed output is m(a, 1) = ^ with an av- 
erage power of 0(a, 1) = 0.5. In this scenario, compression is energy-expensive 
compared to transmission, and, because compression is not required for stabihty, it 
is easy to see from (44) that transmitting aU data without compression is optimal. 
Thus, the poHcy of (44) has k(t) = for all i, and transmits whenever U(t) > 2048, 
yielding an optimal average power Ptot — bE{A{t)} /iimax = 0.5. 

We simulate our dynamic compression and transmission algorithm over 10^ slots, 
for various choices of the V parameter. Fig. 2 shows that average power indeed 
converges to 0.5 as V is increased. Fig. 2 also shows that, as expected, incorporating 
the simple delay improvement of Section 5 does not affect power expenditure. Fig. 
3 shows the time average queue backlog versus V. For simulations without the delay 
improvement, the queue backlog increases linearly with V. The delay improvement 
uses Uthresh = niax[T^//Xmaa; — /imaa; , 0] , and rcduccs queue backlog (maintaining a 
relatively constant average backlog for V > 5000). The dashed horizontal line at 
U — 1920 bits (shown in Fig. 3) is the average queue size obtained by the k{t) 
policy that performs no compression and transmits only when U{t) > fimax- 
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Fig. 4. Avg. Power Consumption vs. Raw Arrival Rate. Data points for the experiments without 
compression are shown only for the stable region, i.e., for raw arrival rates less than 2048. 



7.2 Scenario II: Compression for Stability 

We next consider the same scenario, but increase the raw data rate beyond fimax, so 
that compression is required for queue stabiUty (this removes the third "Singularity 
Assumption"). However, the proper fraction of time to compress may be different 
for each observed A{t) value, and in general it depends on the distribution of the 
arrival process A{t). Our dynamic algorithm optimizes without this statistical 
knowledge, learning the correct actions for each observed A(t) value. 

Fig. 4 shows the increase in average power expenditure for our algorithm (with 
delay improvement) as the raw arrival rate increases. This raw arrival rate is 
increased by adjusting the packet size b from 256 to 1024 (the parameter V is fixed 

to V^ = 10 kbits /unit power, so that Uthresh — T^9Ji[V / ^max — fJ'max,0] ~ 2835 

bits, and the simulation time for each data point is five million slots). Also shown 
is the average power when there is no compression but when the same dynamic 
transmission strategy is used. For arrival rates bK{A{t)} > 1024, compression 
is required for energy efhcicncy, and for b¥,{A{t)} > fimax ~ 2048, compression 
is required for both energy efficiency and stability. For V = 10, our dynamic 
algorithm yields energy efficiency within roughly 0.4% of optimal for the rate region 
tested. For example, when the raw arrival rate is 3400, the optimum is P*^ ~ 1.310 
(achievable by compressing whenever A(t) > 3), and our algorithm achieves Ptot — 
1.314. Because compression reduces data on average by a factor of 2, the maximum 
raw arrival rate that can be stably supported is 2fimax- When blL{A{t)} > 2fimax, 
our algorithm learns to compress all data, leading to an average power expenditure 
of 1.5 (0.5 power units for compression, plus 1 unit for transmission). 

Fig. 5 shows the change in average queue size as the raw data arrival rate is 
increased. Without any data compression, the average queue backlog grows to 
infinity as the raw data rate approaches the vertical asymptote 2048 bits/slot. 
With the data compression option the queue size remains quite flat beyond this 
threshold, increasing at a new vertical asymptote at 4096 bits/slot. 
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Fig. 5. Avg. Queue size (bits) vs. Raw Arrival Rate. 

7.3 Scenario III: Nonlinear Rate-Power Curve 

For this scenario, the Bernoulh arrival process is the same as in Scenario I (Section 
7.1), with packet size b = 256 bits. However, we make the following changes: 

— The rate-power curve is non-linear in power with C{P) — Q!log(l + /3P) for 

transmit power P, where P is any real number in the interval < P < Pmax-^ 
— The raw data arrival rate is less than the maximum transmission rate, i.e. 

bE{A{t)} = 1024 < ^imax=a\0g{l+f3P^ax)■ 

— The compressed data was taken from a trace of experimental data from [Pack 
et al. 2006], and was compressed using the zlib compression library [www ]. 

A single compression option (/C = {0,1}) is available at the transmitting node. 
For fc = 1, we have m{a,k) — j^ for A{t) — a < 3 packets, and m{a,k) — 
j^ for A{t) ~ a > 3. These average compression ratios were obtained from the 
experimental data in [Paek et al. 2006] using the zlib data compression library 
[www ]. The work in [Paek et al. 2006] considers a wireless sensor network where 
each node senses vibrations of a large suspension bridge. We assume the power 
expenditure during compression, (j){a, 1), is 5 units for A{t) — a < 3 and 8 units for 
A{t) = a > 3. For transmission, we use Pmax = 750 power units, a = 1060, and 

(3 = 1/16, so that ^rnax = C {Pmax) ~ 4100 bits/slot. 

Fig. 6 shows that the average power consumption for our dynamic algorithm 
converges to 22.21. As expected, the power curves are almost identical with and 
without delay improvement. Fig. 6 also shows the average power expenditure 
converges to 26.042 if all data is transmitted uncompressed (but still using our 
transmission strategy of (19)). Our dynamic compression algorithm yields an energy 
savings between 15 and 25 percent across the V range tested, as compared to 
sending all the data uncompressed. Fig. 7 shows how the average queue backlog 
increases with V . The plot without delay improvement uses Uthresh = (which 
is also the value of Uthresh that would be given by (36) for the V range tested),^ 



""The log() used here denotes a natural logarithm. 

''The value of Uthresh given by (36) uses l3max = af3 and i^max = C{Pmax)- 
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Fig. 6. Nonlinear Rate-Power Curve: Avg. energy usage vs. V. The two curves with compression 
(with and without delay improvement) are almost identical. 
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Fig. 7. Nonlinear Rate-Power Curve: Avg. Queue size vs. V 

and the plot with delay improvement uses Uthresh as the largest value that satisfies 
Property 1 under the given transmission policy (different from the value given by 
(36)). This value of Uthresh for logarithmic rate-power curves is derived in Appendix 
E (see equation (72)). Note that for a = 1060, Pmax ^ 750, /? = 1/16, V = 70000, 
we have Uthresh = 1056.6 bits, a quite significant improvement in queue backlog 
with no loss of power efficiency. 

8. CONCLUSION 

This paper presents a dynamic decision technique for joint compression and trans- 
mission in a wireless node, using results of stochastic network optimization. The 
approach allows total average power expenditure to be pushed arbitrarily close to 
optimal, with a corresponding delay (and queue congestion) tradeoff. The resulting 
compression and transmission algorithms are simple to implement and operate well 
in a variety of settings. We believe this approach will also be useful for manage- 
ment of compression and sensing in multi-hop networks, where energy, stability, 
and delay issues will become increasingly important in future applications. 
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Appendix A - Proof of Lemma 1 

Proof. (Lcnima 1) The function h*{r) is defined in terms of an infimuni of 
^{Pcomp{t)} over all stationary randomized policies that yield E{R{t)} < r. It 
follows that there exists an infinite sequence of stationary randomized policies, in- 
dexed by integers ie{l,2,...}, having expectations E {i?(*'(t)} and E \ Pcomp{t) \ 
that satisfy: 

E{i?W(t)} <r for alH e {1,2,...} (45) 

lim,_oo E {pi:Lp(t) ]^h* (r) (46) 

However, each policy i is defined in terms of a collection of probabilities {Tak) for 
a S {0,1,..., TV} and fc G /C. This collection of probabilities can be viewed as 
a finite dimensional vector that is contained in a compact set Vt defined by the 
constraints (8)-(9). The compact set VI contains its limit points, and hence the 
infinite sequence {{la k)}?^i contains a convergent subsequence that converges to 
a point (7* ^) G ri. This point is a vector of probabilities that define a stationary 
randomized algorithm with expectations E{i?*(t)} and E|P^*q (i)}. Now recall 
that a general stationary randomized algorithm yields expectations E{i?(t)} and 
E {Pcomp{t)} that can be expressed as linear (and hence continuous) function of the 
probabilities (7a,fe), as shown in the left hand sides of equations (6) and (7). Hence, 
the properties (45) and (46) are preserved in the limit, so that E{i?*(i)} < r and 

If E{i?*(t)} = r, then we are done. Else, we have E{i?*(t)} < r < 6E{A(t)} 
(recall that r < bE{A(t}} by assumption in the statement of Lemma 1). Hence, 
r = eE{R*{t)} + (1 - e)bE{A{t)} for some probability 9. Note that the 0-power 
algorithm of no compression yields an expected compression output of exactly 
bE{A{t)}. It follows that defining R'{t) as the stationary randomized policy that 
chooses R*(t) with probability 9 and chooses no compression with probability {1^9) 
yields E{R'{t)} = r. This policy R'{t) cannot use more power than policy R*{t), 
and hence E{P^'„„p(t)} < h*ir). But we also have h*{r) < E{P^^^p{t)}, be- 
cause h* (r) is defined as the infimum average power over all stationary randomized 
policies that yield a compressor output rate of at most r. D 

Appendix B - Proof of Theorem 1 

Here we prove Theorem 1. Consider any policy that stabilizes the queue, and let k{t) 
and Ptran{t) be the resulting compression and transmission power decisions chosen 
over time (where k{t) G /C and Ptran{t) G V for all t). Let R{t) = ^'{A{t),k{t)) 
be the resulting bit output process from the compressor, and let Pcomp{t) be the 
resulting compression power expenditure process. Let /i(i) — C{Ptran{t)iS{t)) be 
the transmission rate process. We want to show that: 

1 *^^ 

limSUp - ^ E {P,omp{T) + Ptranir)} > P,,*, (47) 

t — ^oo t^ ri 

T=0 

where P*^^ is defined in Theorem 1. We have two preliminary lemmas. 
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Lemma 5. Suppose there are constants r and Pc together with an infinite se- 
quence of times {tij^i such that: 



lim - VE{i?(T)} == r (48) 

i^oo ti ^ — ' 

T = 
1 **"^ _ 

lim - VE{Peomp(T)} - Pc (49) 



t,-.oo t, 

T = 

ThenPc > h*{r). 

Proof. The proof is given in Appendix C. D 

Lemma 6. Suppose there are constants ft and Pt together with an infinite se- 
quence of times {ti}^i such that: 

ti-i 



lim 1 VlE{Pt™„(T)} = Pt (50) 

lim -Y^EMr)} =71 (51) 

* T=0 

Then Pt > g*{ji). 

Proof. The proof is given in Appendix C. D 

Now define Ptot as the lim sup total power expenditure given by the left hand 
side of inequality (47). Let ti be an infinite subsequence of times over which the 
lim sup is achieved, so that: 

1 *--^ _ 

lim - V E {Pcomp(T) + Ptranir)} = Ptot (52) 

t.^oo ti ^^p 

Now define: 

-. t-i -. t-i 

i?(t) = -5IlE{i?(T)} , Pcompit) ^ -Y,^{Pcomp{T)} 

T=0 T=0 

_ t-1 t-1 

Ptranit) = ^ I] ^ {Ptranir)} , W) = '^Y.^ ^^"^^^^ 
r=0 T=0 

and note that for all timcslots t we have: 

< R{t) < bE{A{t)} , < Pcompit) < <j)max 
< Ptranit) < Pmax , < 7l(i) < r^ax 

It follows that iRiti), Pcompiti), Ptraniti),'piti)) Can be viewed as an infinite se- 
quence contained in a four dimensional compact set, and thus has a convergent 
subsequence. Let {ti} represent the convergent subsequence of times, so that there 
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exist constants r, P^^ Pt, and fi such that: 



Um l^E{i?(r)} = 7 
1 *'"^ _ 

Um -y2E{P,ornp{T)}^P, 

1 *'"^ _ 

hm - VE{Pt™„(T)} = P; 

,-»oo t- -"^ ' 

* T=0 

hm -^E{a.(t)} 



T=0 



Furthermore, because {ti}°^i is an infinite subsequence of the original sequence 
{U}, we have from (52) that Pc + Pt = Ptot- From Lemmas 5 and 6 we must have 
the following: 

Pc > h*{r) 

Pt > 9*m 

Therefore: 

Ptot>h*{r)+g*{-p) (53) 

We now use the fact that queue U{t) is stable. It is known that a stable queue 
must satisfy (see [Neely 2006] [Georgiadis et al. 2006]): 

hm M^^ . (54) 

t^oo t 

However, note that for all times ti we have: 

T=0 T=0 

This is true because the total unfinished work in the system at time ti is no more 
than the total bit arrivals minus the maximum possible bit departures during the 
interval from to t^ — 1. Therefore (taking an expectation and dividing by U): 

E{UiU)} ^ *'-' ^ ''-' 



ti 



>l^]E{i?(r)}-l^E{M(r)} 



Taking a limit of the above expression as ti -^ oo, and using (54) yields > r — /I- 
Therefore, queue stability implies that r < JI. Because the function g*{r) is non- 
decreasing, it follows that g*{r) < g*{jl). Using this fact in (53) yields: 

Ptot>h*(r)+g*{r) 

Furthermore, it is not difficult to show that the values of r and /I must satisfy 
Tmin ^ f ^ UE,{A{t)} and < /I < Vmax- Bccause r <JI, it follows that: 

rmin <r< min[rmax,bE {A(t)}] (55) 
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Therefore, the value of h* (r) + g* (r) is greater than or equal to the minimum value 
of this quantity, minimized over all r that satisfy the constraint (55), which is the 
definition of P*^ . Therefore: 

Ptot>h*{r)+g*{r)>P:, 

This proves Theorem 1. 

Appendix C - Proof of Lemmas 5 and 6 

Proof. (Lemma 5) Here we prove Lemma 5. Suppose that (48) and (49) hold. 
For all timeslots t, we have (by iterated expectations): 

E{R{t)} = E{E{i?(t)| A(t),fc(t)}} 
= E{m{A{t),k{t))} 

Similarly, we have E{Pcomp(0} = IE{0(A(t), A:(t))} for all t. The equations (48) 
and (49) thus become: 

ti~i 



lim -V E{m(^(T),fc(T))} = r (56) 

T=0 

lim - VE{0(A(T),fc(T))} - P, (57) 

.- — ^ no r .- ^ — ^ 



ti^oo t, 

' T=0 



For any time t we have: 
t-i 
t 



T = 

N 



i^E{m(A(r),fc(T))} 

N t-1 

J2J2^J2 "^("' k)pA{a)Pr[k{T) = k I A{t) = a] 

a=OkelC T=0 

N 



a=OkelC 

where we define probabilities {"/a.k{t)) as follows: 

t-i 



7a,fc(t)^jE^^[^M^^I^M=«] 



Similarly, for any time t we have: 

t-l AT 

t 



^ t-l AT 

- ^ E mA{T), fc(T))} = 5] 5] 0(a, /c)p^(a)7a,fe(t) 



T=Q a=akeic 
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It follows from (56) and (57) that: 



N 



lim y^y^pA{a)m{a,k)^a.k{'ti) = r (58) 

a=OkelC 

N 

lim y^y^ pA{a)(t>{a, k)ja.k{ti) ^ Pc (59) 

ii^oo ^ — ^ ^ — ^ 
a=Ok<^K 

It is straightforward to show that the probabilities {'^a.k{t)) satisfy the following 
constraints for all t: 

la,kit) > for aUa,fc (60) 

Xl7a,/c(i) = l foraUa (61) 

keic 

The above constraints imply that (7a,fe(i)) can be viewed as a vector of values 
contained in a finite dimensional compact set for all t. It follows that {(7a,fe(ii))} 
forms an infinite sequence of probability vectors contained in a compact set, and 
so there must exist a convergent subsequence of times {t^} for which {(7a,fe(ii))} 
converges to a point (7* j.) contained in the set. Therefore: 

7a,fc > for all a, k (62) 

^7;^ = ! foraUa (63) 

kejc 

N 

^^PA(a)m(a,fc)7*_fc =r (64) 

a=OkGlC 

N 

Yl Yl PAia)^i"'^ ^ha,k = Pc (65) 

a=Ofee/C 

where (62) and (63) follow because (7* j,) is a limit point of the compact set defined 
by (60) and (61) and hence is an element of that set. Equalities (64) and (65) 
follow because {t'^} is an infinite subsequence of the original sequence of times {ii}, 
and hence the same limits in (58) and (59) are preserved when taken over this 
subsequence. 

Because (7* /.) and Pc satisfy (62)- (65), these values define a particular solution 
for the constraints (6)-(9) of the optimization problem of Definition 1 in Section 
3. Therefore, Pc is greater than or equal to the infimum value of power, infimizcd 
over all solutions that satisfy these constraints, which is defined as h*{r). That is, 
Pc > h*{r). This completes the proof of Lemma 5. D 

Proof. (Lemma 6) Here we prove Lemma 6. Suppose that (50) and (51) hold. 
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Similar to the proof of Lemma 5, we can show that for any timeslot t: 



i^EMr)} = ^7r,7I,(t) (66) 

r=0 seS 

1 *"^ - 

-Y,^{Ptran{T)} = ^7r,P,(i) (67) 



T=o ses 

where JI^ (t) and Pg (t) are defined for aU s E S as fohows 

t-i 

t 



7I,(t) 4 jJ2^{CiPtran{T),s)\Sit)^s} 

^ T = 

P,(t) 4 -^E{Pt™„(T)|5(t)-s} 



For each timeslot t and each channel state s G S, the values (/Ig(t), Ps(t)) defined 
above are in the convex hull of the set fig defined below: 

ngA{i^i,p)\per and tx = C{p,.3)} 

The set fig is 2-dimensional. It follows by Caratheodory's theorem [Bcrtsekas et al. 
2003] that any element {jig{t),Ps{t)) contained in the convex hull of ilg can be 
expressed as a convex combination of at most 3 elements of ilg. Thus, there exist 
powers Ps,z{t) e V and probabilities as,z(0 such that: 

3 

{-pgit),Pg{t)) = ^a.,.(t)(C(P,,.(i),5),P«,.(i)) (68) 

2=1 

where X]z=i '^s.z{t) — 1 for all s, t. 

Using (68) and (66), (67), the limit equations of (50) and (51) become: 



ti — >oo 

ses z=i 



lim y^y^TTsag^{ti)C{Pg^z{ti),s) = fi 

3 

lim y^y^ Trsas.z{ti)Ps,z{ti) = Pt 

: — * nn ' ' ' ' 



ti — i-OO ' 

S£S z = l 

It follows that for any e > 0, there exists a stationary randomized policy for 
choosing Ptran{t) as a random function of the observed channel state S{t) such 
that: 

E{CXPtran{t),S{t))} > Jl - 6 
¥.{Ptran{t)] < Pt + e 

This stationary policy can be modified to create another stationary randomized 
policy that has average transmission rate greater than or equal to ]I simply by 
independently choosing Ptran (t) — Pmax every timeslot with some small probability. 
Thus, for any value (5 > 0, there exists a stationary randomized policy for choosing 
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Ptran{t) that yields: 

K{C{Ptranit),S{t))} > JI 

E{Ptran{t)} < Pt + S 

It follows that K{Ptran{t)} In the above policy satisfies K{Ptran{t)} > g*{'p), be- 
cause g*{p) is defined as the smallest average power over the class of stationary 
randomized algorithms that support an average transmission rate of at least JI (see 
Definition 2 in Section 3). Therefore: 

g*m<E{Ptran{t)}<Pt+S 

The above inequality holds for all (5 > 0, and so g*{jl) < Pt, which completes the 
proof of Lemma 6. D 

Appendix D - Proof of Theorem 4 

Here we prove that the Distortion Constrained Compression and Transmission Al- 
gorithm yields performance as given in Theorem 4. Because the algorithm observes 
the current network state and makes control decisions k{t) £ /C, Ptran{t) G V that 
minimize the right hand side of the drift bound given in Lemma 4, we have: 

^{Z(t)) + V¥.{Ptot{t)\Z{t)}<C 

-C/(t)E{C(F;™„(t),^(i))|Z(t)} 
+C/(t)E{m(A(i),r(i))| Z{t)} 

-x(t)EK, -d(A(t),r(i))| z(t)} 

+T/E{P,;„„(i) +0(A(t),r(t)) I Z{t)] (69) 

where k*{t) G /C and Ptrani^) G V are any other feasible control actions for slot 
t. We obtain bounds on U, Ptot, and D using three different k*{t) and Pfrani^) 
policies. 

— (U Analysis): Let Pj*^^„(i) — Pmax, and let fc*(i) be the stationary randomized 
policy that makes decisions independently of the current queue state, and yields 
the minimum output rate rd.min from the compressor, subject to the distortion 
constraint: 

E{m(A(t),r(t))} = rd,„„ 

E{d(A(t),r(t))} < dav 

Plugging this into (69) yields: 

^{Z{t)) + V¥.{Pt,t{t)\ Z{t)} < 

^ C' y^jyniax '^d,niin\ i * [-^jnax i (pmax\ 

and hence: 

A(Z(t)) < C - U{t)[rmax - Tdjnin] + V[Pmax + 4>raax\ 

Using this drift inequality directly in the Lyapunov Drift Lemma (Lemma 2) 
yields the bound on U given in (42). 
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-{P Analysis): Let Ptrani^) ^^^ ^*(0 be the stationary randomized algorithms 
that choose actions independently of queue backlog and yield: 

E{C(P;™„(t),^(t))} = r* 

¥.{m{A{t),k*{t))} = r* 
E{^(A(i),r(i))} = h*,{r*) 
¥.{d{A{t),k*{t))} < da. 

where r* is the optimal solution of problem (39), satisfying: 

h:{r*)+9*{r*) = PL 
Plugging this into (69) yields: 

A{Z{t)) + VE{Ptot{t) I Z{t)} <C + VP:, 

Using this drift inequality in the Lyapunov Drift Lemma (Lemma 2) yields the 
Ptot bound of (41). 

-(D Analysis): Let Pt^ani^) ~ Pmax a-ud let k*{t) be any stationary randomized 
policy that is independent of queue backlog and that yields: 

E{m(A(t),r(t))} < r^ax 
E{d{A{t),k*{t))} ^ dav-e 

for some value e > 0. Plugging this into (69) yields: 

A{Zit)) + VE{Ptot{t)\ Z{t)} < 

C - X{t)e + V[Pmax + 4>max\ 

and hence: 

A(Z(t)) <C- X{t)e + V[P„,ax + 4>max] 

Using this drift inequality directly in the Lyapunov Drift Lemma (Lemma 2) 
yields: 



limsup-^E{X(T)} 

t — *00 L^ 



It follows that the virtual queue X{t) is strongly stable. Because it has a finite 
maximum departure rate davi the time average expected arrival rate to X{t) 
(given by D) is less than or equal to the time average expected transmission rate 
(given by day) [Neely 2006] [Georgiadis et al. 2006]. This proves (43). 

Appendix E - Derivation of Uthresh for the Logarithmic Rate-Power Curve Model 

The logarithmic model has C{P) = a \og{l+(3P) (using a natural log), with fimax — 
alog(l + pPmax)- The dynamic transmission algorithm solves: 

Maximize: U{t)a\og{l + (3P) - VP (70) 

Subject to: < P < P,nax 
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The largest value of Uthresh that satisfies Property 1 is given by Uthresh = max[0, 6] , 
where 9 is the minimum value for the following optimization problem: 

Minimize: 6 = U - ^i*{U) (71) 

Subject to: ^<U<^+ ^f^ 

where ^*(C/) = Q!log(l + fiP*{U)) and P*{U) is the optimum solution to (70) for 
U{t) = [/, given by: 

I ! nj I ^ 

P*{U) 



Ua 1 



where the operator [x]^ is defined [a:;]Q:^max[0, min[x,y]]. This can be understood 
as follows: The value U — IJ*{U) is the queue backlog after transmission when 
U(t) = U. If U{t) < V/{al3) then P{t) = and fi{t) = (so queue backlog 
cannot further decrease) while if U{t) > V/{aP) + VPmax/a then P{t) = Pmax 
and ^{t) = fimax, so the queue cannot drop below the value it would drop to if 
U{t) =. V/{af3) + VPmax/a. 
The problem (71) reduces to: 

Minimize: U — a \og{Ua(3/V) 
Subject to: ^, < U < ^, + ^^^'^^ 

The critical points of the above problem appear at the two endpoints of the interval 
and at the point U — a (ii this point is inside the interval). We thus have: 



Ui 



thr 



max[0, a- a \og{a^l3/V)] if^<a<^ + ^ 



0, min 



al3' al3 ~^ a 



else 



(72) 
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